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ABSTRACT 


his  final  report  covers  the  theoretical  development  of  and  experimental 
results  for  the  stress  analysis  and  strength  characterization  of  thick 
composite  laminates.  Thick  composite  laminates  (consisting  of  38  plies 
or  more)  are  being  evaluated  for  use  In  future  high  performance 
antl-bal 1 istlc  missiles  and  other  aerospace  applications.  The  theoretical 
development  Includes  the  derivation  of  a hlgh-order  theory  of  plate 
deformation  which  accounts  for  the  effects  of  transverse  shear  deformation, 
transverse  normal  strain,  and  a non-linear  distribution  of  the  In-plane 
displacements  with  respect  to  the  thickness  coordinates.  The  theory  Is 
developed  for  both  homogeneous  and  laminated  plates>and  is  presented 
respect ivaly  in  Chapters  1 and  2.  Chanter  \ presents  further  examination 
of  this  high-order  plate  theory  via  stress  solutions  wfTich  have  been  carried 
out  to  assess  its  accurac^F  This  theory.  In  effect,  enables  close 
estimation  of  three-dimensional  stress  components  from  essentially 
two-dimensional  analysis. 

The  effects  of  the  three-dimensional  state  of  stress  on  the  strength  of  a 
thick  laminates  are  analysed  with  a 3-0  failure  criterion.  trt'Chapter  A, 
j£he  tensor  polynomial  method  is  extended  with  full  account  given  to 
three-dimensional  stress  state  effects  and  is  presented  together  with 
experiment  and  evaluation  of  the  coeff icients. combined  theoretical 
development  of  thick-plate  analysis  and  three-dimensional  failure  criterion 
is  expected  to  improve  confidence  and  full  utilization  of  composites  in 
applications  where  the  thickness  of  the  laminates  prevents  them  from  being 
adequately  treated  by  current  thin-plate  formulations. 
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SUMMARY 


This  final  report  covers  the  theoretical  development  of  and 
experimental  results  for  the  stress  analysis  and  strength  characterization 
of  thick  composite  laminates.  Thick  composite  laminates  (consisting 
of  38  plies  or  more)  are  being  evaluated  for  use  In  future  high 
performance  ant  I -bal 1 1 s t Ic  missiles  and  other  aerospace  applications. 

The  theoretical  development  Includes  the  derivation  of  a high-order 
theory  of  plate  deformation  which  accounts  for  the  effects  of  transverse 
shear  deformation,  transverse  normal  strain,  and  a non-linear 
distribution  of  the  in-plane  displacements  with  respect  to  the 
thickness  coordinates.  The  theory  is  developed  for  both  homogeneous 
and  laminated  plates  and  is  presented  respectively  In  Chapters  I and 
2.  Chapter  3 presents  further  examination  of  this  high-order  plate 
theory  via  stress  solutions  which  have  been  carried  out  to  assess 
Its  accuracy.  This  theory,  in  effect,  enables  close  estimation  of 
three-dimensional  stress  components  from  essentially  two-dimensional 
analyst  s. 

The  effects  of  the  three-dimensional  state  of  stress  on  the 
strength  of  a thick  laminates  are  analysed  with  a 3'0  failure 
criterion.  In  Chapter  k,  the  tensor  polynomial  method  Is  extended 
with  full  account  given  to  three-dimensional  stress  state  effects  and 
Is  presented  together  with  experiment  and  evaluation  of  the 
coefficients.  The  combined  theoretical  development  of  thick-plate 
analysis  and  three-dimensional  failure  criterion  Is  expected  to 
Improve  confidence  and  full  utilization  of  composites  in  applications 
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wher«  the  thickness  of  the  lamlnetes  prevents  them  from  being  adequately 
treated  by  current  thln-plate  formulations. 


In  many  structural  applications,  special  characterization  methods 
are  needed  to  model  very  thick,  laminated  composites.  These  thick 
composites  are  currently  being  evaluated  for  future  high  performance 
antibal I istic  missiles  and  other  aerospace  applications.  A rational 
characterization  must  include  1)  extending  current  two-dimensional 
laminated-plate  theory  to  reflect  the  three-dimensional  characteristic 
of  these  thick  plates  for  the  purpose  of  accurately  estimating 
deflection  and  the  three-dimensional  state  of  stress,  2)  establishing 
failure  criterion  to  assess  the  effect  of  the  three-dimensional  state 
of  stress  and  hence  of  the  strength  of  a thick  laminate.  In  light  of 
these  requirements,  a high-order  plate  deformation  theory  is  derived 
to  remove  the  restriction  of  the  plane-sectlon-remains-plane 
hypothesis  in  conventional  laminated  plate  theory  by  allowing  a 
plane  section  to  assume  third-order  displacement  modes.  This  theory 
is  fully  developed  for  homogeneous  and  laminated  plates  and  is 
presented,  together  with  stress  solutions,  in  chapters  1,  2 and  3- 
In  order  to  present  the  theory  in  complete  and  concise  form  the 
theories  and  stress  solution  are  presented  in  the  format  of  self- 
contained  chapters.  Finally  the  effects  of  a three-dimensional  state 
of  stress  on  the  strength  of  the  thick  laminates  are  assessed  with 
a three-dimensional  tensor  polynomial  failure  criterion  in  chapter  A. 
These  four  chapters  together  provide  the  necessary  tools  for  the 
stress  analysis  and  strength  characterization  of  thick  laminates 
composites.  Thesn  tools  should  enhance  a more  complete  and  confident 
utilization  of  composites  in  the  form  of  thick  laminates. 


Chapter  1 


A HIGH  ORDER  THEORY  OF  PLATE  DEFORMATION 
Homogeneous  Plates 


K.  H.  Lo,  R.  M.  Christensen,  and  E.  M.  Wu 


ABSTRACT 


A theory  of  plate  deformation  is  derived  which  accounts  for  the 
effects  of  transverse  shear  deformation,  transverse  normal  strain, 
and  a nonlinear  di str ibution  of  the  in-plane  displacements  with  respect 
to  the  thickness  coordinate.  Tha  theory  is  compared  with  lower  order 
plate  theories  through  application  to  a particular  problem  involving 
a plate  acted  upon  by  a sinusoidal  surface  pressure.  Comparison  Is 
also  made  with  the  exact  elasticity  solution  of  this  problem.  It  is 
found  that  when  the  ratio  of  the  characteristic  length  of  the  load 
pattern  to  the  plate  thickness  is  of  the  order  of  unity,  lower  order 
theories  are  inadeouate  and  the  present  high  order  theory  is  required 
to  give  meaningful  results.  The  present  work  treats  homogeneous 
plates  while  Chapter  II  involves  laminated  plates. 

INTRODUCTION 

The  development  and  application  of  classical  plate  theory  is  one 
of  the  achievements  of  modern  engineering.  It  is  continuously  being 
applied  to  new  problems  to  gain  new  and  needed  design  information. 

Despite  its  successes,  however,  the  inherent  limitations  of  the  classical 
theory  necessitate  the  development  of  more  refined  and  higher  order 
theories  of  plate  behavior.  More  sophisticated  models  of  plate  behavior 
find  application  to  problems  where  classical  plate  theory  is  simply 
inadequate  to  describe  the  behavior.  Such  examples  concern  plate  with 
cutouts,  contact  problems  involving  plates,  and  laminated  plates.  The 
present  work  concerns  the  derivation  and  evaluation  of  a particular  high 
order  theory  of  plate  behavior. 
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Before  describing  the  present  theory,  it  is  necessary  to  briefly 
review  the  recent  developments  in  the  generalization  of  classical 
plate  theory.  Reissner  [1,2]  was  the  first  to  provide  a consistent 
theory  which  incorporates  ♦he  effect  of  shear  deformation.  The  derivation 
given  by  Reissner  resulting  in  displacements  of  the  form 

u ■ u#  + zto 
rx 

v - v8  + z<J/y  (1) 

w ■ w° 

where  z is  the  coordinate  normal  to  the  middle  plane,  and  u° , v°, 

V , V , and  w°  have  a dependence  upon  the  in-plane  coordinates  x and  y, 

x y 

and  , ip  , and  w#  are  actually  weighted  averages.  The  basic 
x y 

assumption  used  by  Reissner  involved  consistent  forms  for  the  stress 
distributions  across  the  thickness.  A special  variational  theorem  was 
used  to  determine  both  the  equations  of  equilibrium  in  terms  of 
resultants  and  the  stress-strain  relations  in  the  form  involving 
resultants  and  the  functions  in  (1).  At  the  same  level  of  approximation, 
Hindi  in  [3]  employed  kinematic  assumptions  of  the  form  of  (1),  and 
without  introducing  corresponding  stress  distribution  assumptions, 
obtained  the  governing  equations  from  a direct  method.  In  Hindi  in's 
derivation  it  was  necessary  to  introduce  a correction  factor  into 
the  shear  stress  resultants  to  account  for  the  fact  that  relations  (1) 
predict  a uniform  shear  stress  through  the  thickness  of  the  plate, 
which  is  incorrect  and  in  general  would  violate  surface  conditions. 

The  correction  factor  was  evaluated  by  comparison  with  an  exact 
elasticity  solution.  It  is  useful  to  observe  that  the  form  (1)  applies 
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to  both  the  classical  theory  of  plate  bending  as  well  as  to  the 
theories  of  Reissner  and  Mindlin  which  include  the  effect  of  transverse 
shear  deformation.  Thus,  considering  the  terms  in  (1)  as  the  first 
terms  in  a power  series  expansion  in  z,  it  is  seen  that  the  classical 
theory  and  the  shear  deformation  theory  are  of  the  same  order  of 
approximation.  The  classical  theory  is  merely  a special  case  of  the 
shear  deformation  theory,  wherein  the  shear  modulus  in  terms  associated 
with  the  transverse  shear  deformation  is  taken  to  be  very  large, 
such  that  transverse  shear  deformation  can  be  neglected. 

There  have  been  several  theories  proposed  which  are  of  higher  order 
than  those  of  Reissner  and  Mindlin.  Typical  efforts  along  these  lines 
will  be  mentioned  here.  The  next  higher  order  theory  from  that 
embodied  in  (1)  involves  displacement  forms  of  type 

u ■ u°  + zip 

x 

v - v°  + zip  (2) 

2 

w * w°  + zil)  + z C 
z z 

which  invludes  the  effect  of  transverse  normal  strain.  Displacement 
assumptions  of  the  form  of  (2)  along  with  corresponding  stress 
distribution  assumptions  have  been  used  by  Naghdi  [4]  to  derive  a 
general  theory  of  shells,  and  by  Essenburg  [5]  to  derive  the  corresponding 
one-dimensional  plate  theory.  In  the  context  of  contact  problems, 

Essenburg  [5]  demonstrated  the  utility  and  advantages  of  the  theory 

II 

based  upon  (2)  over  lower  order  theories.  Whitney  and  Sun  [6]  also 
utilized  assumptions  of  the  level  of  (2)  to  develop  a theory  of 
lamianted  cylindrical  shells.  However,  there  is  an  inconsistency  in 
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their  approach.  They  used  a shear  correction  factor  of  the  same  type 
as  that  employed  by  Hindi  in  in  deriving  stress  resultants.  Whereas 
a factor  of  this  type  was  appropriate  to  Hindi  In's  derivation  since  It 
assumed  uniform  shear  stresses  across  the  thickness,  the  same  type  of 
correction  factor  is  not  appropriate  for  use  with  the  displacements  of 
the  form  of  (2).  This  Is  because  non-uniform  shear  stresses  are 
implied  by  (2)  along  with  consequent  satisfaction  of  top  and  bottom 
boundary  conditions  of  shear  tractions;  thus  the  rationale  for  a 
correction  factor  is  obviated. 

The  next  higher  level  theory  is  based  upon  the  assumed  displacement 

forms 

u - u°  + + z2;x 

v - v°  + z<|i  + (3) 

2 

w - w®  + z4»z  ♦ z 

A theory  derived  from  (3)  has  been  given  by  Nelson  and  torch  [7]  for 
application  to  laminates.  This  theory  however  has  the  same  defect 
in  application  as  that  mentioned  above  In  connection  with  Ref.  [6]; 
namely,  a shear  correction  factor  was  employed  when  in  fact  it  Is 
inconsistent  with  the  level  of  approximation  in  (3).  Hildebrand, 
Reissner  and  Thomas  [8]  briefly  examined  a theory  of  the  level  of 
(3)  and  concluded  that  the  inclusion  of  the  quadratratlc  terms  in  the 
in-plane  displacements  dees  not  provide  a significant  advantage  over 
the  lower  level  theory. 

Reissner  [9]  has  presented  a theory  which  to  a consistent  degree 
of  approximation  gives 
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U • z<l>  + Z3<J> 

X X 

V - z(J)y  + z\  (4) 

w ■ w#  + z^C 

z 


where  the  last  relation  in  (4)  follows  by  combining  relations  (9)  and 
(11)  of  Ref.  [9].  Reissner  has  shown  that  the  plate  theory  corresponding 


to  (4)  gives  very  accurate  results  compared  with  the  elasticity  solution 
for  the  bending  of  a plate  with  a circular  hole.  It  should  be  noted 
that  the  theory  based  upon  (4)  represents  the  lowest  order  correct /on 
for  out-of-plane  deformation  effects  to  the  classical  theory  embodied 
in  the  first  terms  in  (4).  Though  these  results  obtained  by  Reissner 
are  impressive,  a theory  based  upon  (4)  neglects  the  contribution  of 
in-plane  modes  of  deformation;  only  out-of-plane  effects  are  considered. 
Such  in-plane  effects  may  be  of  importance  in  certain  plate  problems, 
and  this  effect  will  be  investigated  herein. 

The  theory  to  be  presented  here  is  appropriate  to  the  following 
displacement  forms: 


u ■ u°  + zip  + z2Z  + z34> 

X X X 

V - v°  + zip  + Z2z  + z34> 
y y y 

2 

w * w°  + z<J>z  + z Cz 


(5) 


which  is  of  the  same  level  as  the  Reissner  theory  corresponding  to 
(4)  but  includes  both  in-plane  and  out-of-plane  modes  of  deformation. 
The  theory  of  plate  behavior  based  upon  (5)  will  be  derived  by 
application  of  the  principle  of  stationary  potential  energy.  The 
accuracy  of  this  theory  will  be  assessed  by  direct  comparison  with  an 
exact  solution  from  the  theory  of  elasticity. 
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Before  turning  to  the  derivation,  it  is  pertinent  to  outline  the 
motivation  for  the  present  work.  The  primary  intended  application  for 


the  present  high  order  plate  theory  is  in  the  field  of  laminates. 

It  is  well  known  that  laminated  plate  behavior  provides  a particularly 
critical  test  of  the.  Bernoulli  hypothesis  concerning  plane  sections 
remaining  plane.  The  mismatch  in  properties  causes  deviations  from  the 
lowest  order  terms  in  displacement  forms  (4).  However,  to  place  the 
present  work  in  its  proper  context,  it  is  useful  to  derive  it  first  in 
the  form  directly  suitable  for  application  to  homogeneous  plates. 
Accordingly,  the  present  paper,  Chapter  I is  concerned  with  homogeneous 
plate  behavior  while  Chapter  II,  following,  is  concerned  with  the 
application  to  laminates.  The  present  work  thus  affords  the  opportunity 
to  assess  the  importance  of  the  full  form  of  displacements  (5)  compared 
with  the  partial  form  (4)  for  certain  types  of  problems.  This  same 
question  will  be  further  explored  in  the  context  of  laminates  in 
Chapter  I I . 

Finally,  it  should  be  mentioned  that  a theory  of  the  level  of  (5) 
certainly  is  of  a rather  complicated  form,  and  the  question  arises 
whether  there  is  a practical  need  for  such  a theory.  It  is  a question 
of  the  degree  of  accuracy  required.  For  problems  which  involve 
rapidly  fluctuating  loads  with  a characteristic  length  of  the  order 
of  the  thickness,  the  results  presented  herein  show  that  a theory  of 
the  degree  of  sophistication  of  (5)  is  required  to  give  meaningful 
results.  Also,  It  should  be  recognized  that  the  theory  to  be  developed 
here  is  amenable  to  numerical  integration  with  respect  to  the  planar 
x,y  coordinates. 
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THEORY 


Plate  theories  can  be  developed  by  expanding  the  displacements 
in  power  series  of  the  coordinate  normal  to  the  middle  plane.  In 
principle,  theories  developed  by  this  means  can  be  made  as  accurate 
as  desired  simply  by  including  a sufficient  number  of  terms.  In 
practice,  however,  a point  of  diminishing  returns  is  reached  whereby 
the  complexity  of  the  resulting  forms  becomes  too  great.  Here,  we 
seek  the  minimum  number  of  terms  which  include  the  effects  of 
transverse  shear  deformation,  transverse  normal  strain,  and  warpage 
of  the  cross  section.  Thus  the  displacements  are  taken  in  the  form 
of  (5). 

The  principle  of  stationary  potential  energy  is  used  to  derive 
the  governing  equilibrium  equations.  It  is  found  that 


N + N 

+ q “0 

x,x  xy,y 

X 

N + N 

+ q *0 

y,Y  xy.x 

Y 

Q + Q -* 

x,x  >,y 

• q = 0 

M + M 

- Q + m - 0 

x,x  xy,y 

X X 

M + M 

y,y  xy,x 

" ^Y  + mY  " ° 

R + R 

• N + m - 0 

x ,x  y,y 

z 

p + P 

- 2R  + n = 0 

x , x xy,y 

X X 

p + P 

- 2R  + nw  - 0 

Y,y  xy.x 

Y Y 

S + S 

- 2H  + n ■ 0 

x , x y,y 

z 

M + M „ 

“ 3S„  + * ■ 0 

x,x  xy.y 

X X 

M + M 

" 3S  + 1 - 0 

Y.y  xy.x 

y y 

(6) 
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where  the  stress  resultants  are  defined  by 


N N N N Q Q 
x y i xy  * Y 


M M H M R 
x y z xy  x y 


p p p 
X y xy 


Mx  My  Mxy 


h/2  P 


-h/2  l z 


tox  °y  °z  Txy  Txz  Tyz]dz  (7) 


h/2  r * 


-h/2  L z 


(o  a t ) dz 
x y xy 


rs  S 1 - I z2  ( T r ) dz 

iav  SyJ  xz  yz 

y j -h/2 


(qx  nx)  - [Tx2(h/2)-txz(-h/2)l  [1  h A] 

(qy  ny)  - CTy#(h/2)-Tyl(-h/2J]  [I 
(mx  lx)  - [txz(h/2)  +rxz(-h/2)1  [h/2  h3/8] 

(my  iy)  - [Tyz(h/2)  + Tyz(-h/2)l  [h/2  h3/8] 

(q  n)  -[az(h/2)-  oz(-h/2)  [1  *2/4] 

m - \ taz(h/2)  ♦ az(-h/2)l 


The  resultants  in  terms  of  the  displacement  functions  are  given  by 
Nx  - (X+2u)hu* , x + Xhv* ,y  + Xh<|*z  + yj  [(X+2uKx>x  + Xi*y»y^ 

N - Xhu* , + (X+2p)hv* , + Xh4»z  + JJ  [XCx,x  * {U 

y X ' 

h3 

N - Xhu*.  + Xhv*  + (X+2u) hij»z  ♦ jy  I^,x+XCy.Y] 


„Mu\  *V.xl*%UXiY*CyiJtl 


Mx  = T2  [(X+2u)^,x  + A^y,y  + 2AS} 


+ 80  t(A+2u)<l>x,x  + Hy,y] 


My  “ IT  [A*x,x  + {7'+2y)«y.y  + 2A?z] 


+ 80  I*x.x  + (A+2u)\,y] 


Mz  " 71  t^x.x  + A*y,y  + 2<A+2^2] 


+ 80  tA*x,x  + Hy,y] 

Mxy  " % ^x,y  + Vx1  + % lK,y  + ^y,x] 


p = TT  C (X+2y)  u*  , + Xv®  , + \ip  ] 

x 12  x ’y  rz 


+ "So  l(x+2»)5„  „ + J 

3 


x , x ’y,y‘ 

P.  -tt  [Xu*,v  + (X+2u)v®,„  + X<pj 


h 

y 12 
.5 


+ M ^A?-  - + ^A+2y^'< 


x,x 


’y.y 


\ ■ So  »X+»>*X.X  * XVy  + “'z1 


+ m [{X+2uUx,x  + ^y.y1 
\ ■ S5-  [A^x.x  + <A+2^y,y  + 2A*z] 

+ W [A*x,X  + (A+2^^y,yl 

Mxy  " 55”  [*x,y  + Vx1  + ^x,y  + Vx1 
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Qx-UhWx  t3«x»tI>x] 

uh3 

Qy  - yh[^y  + W°,y]  + -{y  [34>y  + CZfyl 


yh 


R -J!h  [2c  1 

x 12  sx  rz,xJ 


12 


[2<y  * *z.y] 


and 


(15) 


(16) 


yh3 

12 

w\x] 

+ uh5 
*~5o 

t3«x 

* 'z.*1 

12 

[*y  * “*.yl 

80 

[3*y 

* *z.y] 

(17) 

where  X and  u are  the  Lame  constants. 

Finally,  the  boundary  conditions  along  the  edge  of  the  plate 
require  that  one  member  of  each  of  the  following  eleven  products  must 
be  presecribed: 


"V  V “V  "n 
V V "n  V %t  ?t 
"nt  V V 

and  S„  C2 

where  n and  t are  the  directions  normal  and  tangential  to  the  edge 
of  the  plate. 
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An  approximate  theory  can  be  critically  assessed  in  comparison 
with  an  exact  result.  Fortunately,  there  are  exact  solutions  available 
from  the  theory  of  elasticity  which  are  suitable  for  the  present 
purposes.  The  solution  of  use  here  is  that  of  the  deformation  of  an 
infinite  plate  of  thickness  h subjected  to  a pressure  on  the  top 
surface  z ■ h/2  of  the  form 

q * qo  sin  T-  08) 

with  all  other  surface  tractions  vanishing  identically.  From 
equations  (1)  it  follows  that 

qx  “ qy  * mx  ” my  " nx  ” ny  * ^x  * % * 0 09) 

with 

h , TTX 

m “ i qo sfn  *r 
2 

h • TTX  f 

n * ip  qQ  Sin  — (20) 

The  equations  of  equilibrium  (6)  take  the  following  special  forms 
for  this  problem 
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’'■‘"’x.x  * * % '5*x,x  * C,,»]  *1-0 

T7  m+2u)*x,xx  * 2X'z.x>  * SS  [(A*2u)*x  „] 

-Kht**  ♦ *\„)  - ^ (3*x  * CI  X)  ■ 0 
h3 

77  C(u-2X)*  + yw*  - 4(X+2u)C  ] 


* So  I(J“-2X>*x,x  * «z,xx'  * " ■ ° 


Bo  xx  + (2A”3u)c,  _ - 9u$J 


z , x x 


-ft  rr\j.T..\,k  i _ yh3 


+ Utf  t(x>2u)«j.XtXXl  5-  N»x  ♦ w\  1 - 0 


(X+2y)hu° >xx  ♦ ^ §x  + (X+2y)  jj  Cx>xx  - 0 


— I(2u-X)r,^  x «■  U'J'Z  XX1  - Xhu’*x  " (x+2y)h<i»z  + m - 0 


IT  I(x+2w)u\x  ' *«x  + ^^z.x1 


+ S0  I(A+2u)<x.xx>  « 0 


yhv* , + C “ 0 

xx  )2  y,xx 

h3  . , h3  . . uh5 


TT  v*’xx  ' <l‘1  TT  S tl!So  'y,xx  * 0 


and 


L 


4 *»„  * 4 V~  ■ uh*y  - 4 


$ - 0 
y 


•5—  [yij/  -gy$  ] + - ii£L_  m 0 

80  y,XX  J^yJ  1*48  vy,xx  4 Vy  u 


Uh' 


Note  that  the  equations  partially  decouple  such  that  the  equation 
grouping  in  each  of  the  above  number  sets  are  individually  determinate. 

Obviously,  the  solution  for  the  generalized  displacement  functions 
involve  terms  proportional  to  sin  irx/l  and  cos  ttx/L.  Making  the 
appropriate  assumptions  of  this  type,  relations  (21)  - (24)  can  be 
solved  analytically  in  closed  form. 

The  completed  solution  can  be  used  to  evaluate  the  displacements 
and  stress  for  comparison  purposes.  The  mid-plane  displacement  is 
given  by 


w*  - ip  33600  (l-v)(l-2v)  + [7200-16920V  + 5520v2]  tt2 

+ IAOvO-vJit1* 


o-v)V(£)6 


s 1 n 


7TX 


[' 


2 •% 

4( l-v)  | 8400  (l-2v)  + 120  (l-v)ir2  + (l-v)V  | 


17 


(24) 


(25) 


# 


J 


The  solution  for  the  stress  component  a*  is  given  by 


15  (1-v)2tt2  (h/L)2  [12  (z/h)2-l 


* ° [(720  (l-2v)  + 24  (1-v)*2  (h/L)2  + (l-v)V  (h/L)M 

^ 2 2 2 

+ [^2  (k)  ^ f3360°  (,"2V)  * 120  (,'v)  0°-7v) 


- 80  (1-v) 


+ 12(l-v)  |j=-j  3 [2300  (2-v)  + 280  (I-v)tt2  2 * ] | sin  J5- 
j [16800  (l-2v)  + 240  (1-v)  ir2  ^ 2 (1-v) V k ] 


These  results  are  to  be  compared  with  the  exact  solution,  taken  from 
Little  [10],  Other  lower  order  approximate  theories  also  will  be  compared 
with  the  present  results.  First,  the  classical  theory  result  is  noted  to  be 
.4  /_  \ 


"•-T  V •"(¥) 


The  shear  deformation  Reissner  plate  theory  produces  the  result 


• qo  L*1  T,  u2  (2-v)  / h \2  ] nx 

w mr  -T  ['  +To  • tT^rrr]  J sm  - 


Essenburg's  theory  gives 
<J  . *» 


w • - -o  k-  r,t/.  t.2-v)  ( h Y , * . v ( h \ 

D L 10  (l-v)  U J 10  TRT  [ L J 


/ h V*  ( h \ ^ 1 . ttx 

+ To  irj  ■ mo  [l]  Js,n  T 


(29) 


Finally,  for  comparison,  the  present  theory  with  £ m Z - $ »(J>  -0 

x y X y 


in  (5)  has  been  evaluated  and  will  be  referred  to  as  the  "Level  (2)  Theory", 
consistent  with  equations  (2).  This  solution  is 


• flo  l4  r, 

“ T JT  L'  ♦ 


20tt 


(1-v)  [120  + (I-v)tt  ( h/L) 2 ] 


(*)’ 


(30) 


(31) 


DISCUSSION 


Through  comparison  with  an  exact  solution,  we  seek  to  determine 
the  relative  accuracy  of  the  various  approximate  theories,  including  the 
one  presented  herein.  Due  to  the  high  order  of  the  terms  included  in 
the  present  theory,  it  is  of  course  not  convenient  to  use.  Accordingly, 
it  would  be  helpful  to  deduce  guidelines  by  which  on  can  ascertain 
when  it  is  necessary  to  use  a high  order  theory,  as  given  here,  and 
when  a lower  order  theory  will  suffice. 

The  range  of  theories  to  be  considered  are  from  the  classical 
case  to  the  present  form.  In  between  these  extremes,  the  well-known 
Reissner  shear  deformation  theory  will  be  considered  along  with  the 
extension  of  i t by  Essenburg  to  include  transverse  normal  strain 
effects.  Essenburg 's  theory  employs  assumptions  upon  both  stresses 
and  displacements.  The  Level  (2)  theory  referred  to  in  the  previous 
section  uses  exactly  the  same  assumuptions  upon  displacement  as  in 
Essenburg' s theory;  however,  no  corresponding  assumptions  are  made 
upon  the  stresses.  Rather,  the  Level  (2)  theory  solution  is  obtained 
as  a special  case  of  that  given  herein,  which  of  course  is  derived 
directly  from  potential  energy.  Thus,  the  Essenburg  theory  and  the 
Level  (2)  theory  are  of  the  same  order,  but  involve  different 
derivations,  and  it  will  be  of  interest  to  compare  them  directly  with 
each  other. 

Stress  distributions  across  the  thickness  of  the  plate  are 
displayed  in  Figure  1 for  h/L  » 1.5*  where  L is  the  half  wave  length 
of  the  sinusoidal  loading  pattern.  Thus,  the  ordinate,  h/L,  is  the 
ratio  of  thickness  of  plate  to  characteristic  length  of  the  loading 
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pattern.  The  high  order  theory  due  to  Reissner  corresponding  to 
Level  (4)  is  of  the  same  order  as  the  present  theory  based  upon  (5), 
the  difference  being  that  Reissner's  theory  omits  the  combined  effects 
between  in-plane  and  out-of-plane  deformation  modes.  Thus,  the  Reissner 
Level  (4)  theory  of  necessity  produces  antisymmetric  stress  states 
about  the  middle  plane  of  the  plate.  From  Fig.  1 It  is  seen  that  the 
exact  solution  deviates  strongly  from  the  ant i -symmetry  characteristics 
just  mentioned;  thus  it  is  clear  that  a theory  of  the  type  of  (4) 
could  not  come  close  to  reproducing  the  exact  solution.  Also,  it  is 
obvious  from  Figure  1 that  the  Reissner  Level  (1)  and  the  Essenburg 
theories  provide  no  improvement  over  the  classical  theory  in  terms 
of  the  accuracy  of  the  stress  representation.  Considering  the  complex 
shape  of  the  exact  solution  stress  distribution,  the  present  high  order 
theory  provides  an  effective  modeling  result. 

The  maximum  value  of  stress  a (flexural  stress)  Is  plotted 
against  h/L  in  Figure  2.  Clearly,  when  the  ratio  of  the  thickness  to 
the  characteristic  length  of  the  load  is  of  the  order  of  1,  the  present 
high  order  theory  is  needed  to  properly  model  the  deformation  effects 
in  the  plate  and  lower  order  theories  are  inadequate.  This  specific 
conclusion  of  course  applies  only  to  the  present  example,  but  we 
speculate  that  in  all  problems  where  disturbing  features  have  a 
characteristic  length  of  the  order  of  the  thickness,  then  a theory  at 
least  of  the  order  of  the  present  one  would  be  required  to  properly 
model  effects. 

The  maximum  displacement  of  the  middle  plane  of  the  plate  according 
to  the  various  theories  are  compared  with  the  exact  result  in  Figure  3. 
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It  is  seen  that  for  h/L  - 1.5  the  deviations  between  the  approximate 
theories  and  the  exact  result  are  substantial,  and  the  deviations 
continue  to  increase  with  increasing  h/L.  The  present  high  order 
theory  is  seen  giving  the  result  closest  to  the  exact  result.  The 
results  from  Essenburg's  theory  and  the  Level  (2)  theory  are  of  about 
the  same  level  of  accuracy  as  one  would  expect  for  the  reasons  described 
above.  Note  also  that  the  present  high  order  theory  is  clearly  superior 
to  the  Level  (2)  theory. 

The  main  result  of  the  present  work  is  viewed  as  a complement  to 
the  result  found  by  Reissner  [9],  this  result  being  that  a theory  of 
the  level  of  (4)  or  (5)  is  needed.  In  general,  to  provide  a significant 
improvement  in  the  level  of  accuracy  over  that  afforded  by  the  classical 
theory  of  plate  deformation.  Further,  it  follows  that  with  a theory 
of  the  level  of  (4)  or  (5),  very  accurate  results  can  be  obtained. 

It  is  now  possible  to  answer  the  question  raised  in  the  Introduction 
of  whether  in  a given  problem  the  coupling  effects  implicit  in  (5)  are 
needed  rather  than  using  the  simpler  form  (4).  In  the  example  studied 
by  Reissner  [9]  of  the  bending  of  an  infinite  plate  with  a circular 
hole,  the  results  derived  from  a theory  corresponding  to  (4)  were 
entirely  satisfactory  and  sufficient,  whereas  in  the  present  example 
the  theory  based  upon  (5)  is  required  for  the  reasons  described  above. 

An  examination  of  the  governing  set  of  differential  equations  (6)  shows 
that  the  sets  of  equations  governing  the  In-plane  and  the  out-of-plane 
deformation  modes  completely  decouple.  In  the  problem  studied  by 
Reissner  [9],  the  in-plane  contributions  to  the  problem  are  easily 
shown  to  vanish  Identically,  thus  a theory  based  upon  (5)  provides  no 
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new  information  over  one  based  upon  (A)  In  that  problem.  However,  in 
the  present  problem  the  in-plane  contribution  to  the  solution  has  been 
shown  to  be  significant,  and  cannot  be  neglected.  Thus,  it  is  seen 
that  for  plate  bending  problems  where  the  loading  characteristics  possess 
a high  degree  of  asymmetry,  with  respect  to  the  middle  plane,  then  a 
theory  of  the  type  of  (5)  rather  (A ) is  required,  while  problems  with 
loading  characteristics  which  are  close  to  being  anti -symmetrica I with 
respect  to  the  middle  plane  then  a theory  of  the  type  of  (A)  is 
sufficient.  The  bending  of  an  infinite  plate  with  a circular  hole 
is  an  example  of  the  latter  type  of  problem,  while  the  problem 
considered  herein  as  well  as  contact  problems,  are  examples  of  the 
former  type.  Of  course,  neither  type  of  high  order  theory  is  needed 
unless  the  disturbing  feature  of  geometry  or  the  characteristic  load 
length  are  of  the  order  of  the  plate  thickness. 

Finally,  it  is  appropriate  to  mention  the  nature  of  the  stress 
resultants  involved  in  the  present  theory.  As  seen  from  (7)  - (9), 
there  are  resultants  of  higher  order  than  bending  moments  and  shear 
force  involved  in  the  governing  differential  equations.  These  same 
higher  order  resultants  are  also  necessarily  involved  in  the 
specification  of  edge  conditions.  It  is  logical  to  ask  what  is  the 
significance  of  such  high  order  resultants  in  edge  condition 
specifications,  and  is  there  any  way  to  avoid  involvement  with  them. 

The  answer  to  this  question  is  very  simple.  There  is  no  way  to  avoid 
consideration  of  these  high  order  resultants  in  the  present  context; 
indeed,  it  would  be  disturbing  if  the  present  high  order  theory  did 
not  require  the  specification  of  corresponding  high  order  resultants 
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along  edges.  For  traction-free  edges,  of  course,  the  resultants  of  all 
orders  simply  vanish.  For  loaded  edges,  the  distribution  of  tractions 
across  the  thickness  must  be  obtained  or  assumed,  from  which  the 
resultants  of  all  orders  can  be  determined. 

The  extension  of  the  present  plate  theory  to  laminated  plate 
conditions  Is  of  particular  Importance  since  It  Is  known  that  for 
laminates  the  distribution  of  in-plane  displacements  across  the  thickness 
may  be  strongly  nonlinear.  This  extension  Is  presented  in  Chapter  II  of 
this  same  work. 
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Chapter  2 


A HIGH  ORDER  THEORY  OF  PLATE  DEFORMATION 
Laminated  Plates 


K.  H.  Lo,  R.  M.  Christensen,  and  E.  M.  Wu 


ABSTRACT 


The  high  order  theory  of  plate  deformation  developed  in  Chapter  I of 
this  work  is  extended  here  to  model  the  behavior  of  laminated  plates. 
Through  comparison  with  elasticity  solutions,  it  is  shown  the  present 
theory  correctly  models  effects  not  attainable  from  the  classical  theory. 


INTRODUCTION 

With  the  increasing  use  of  composite  materials  in  thick  laminated 
form,  the  need  for  advanced  methods  of  analysis  is  obvious.  For  such 
laminated  systems,  the  components  of  stress  and  strain  transverse  to 
the  plane  of  the  laminate  strongly  influence  the  behavior.  Thus, 
classical  laminated  plate  theory,  which  is  not  formulated  to  account 
for  the  effect  of  these  transverse  stress  and  strain  components,  is  not 
applicable  to  thick  laminates.  A high  order  theory  of  plate  behavior 
is  herein  developed  for  application  to  laminates;  this  theory  is  an 
extension  of  that  developed  in  Part  I of  this  work  [1],  for  application 
to  homogeneous  plates. 

Many  different  high  order  laminated  plate  theories  have  been 
proposed  which  are  intended  to  improve  upon  the  classical  laminated 
plate  theory  by  accounting  for  the  effects  of  the  transverse  components 
of  strain  in  the  plate.  Typical  examples  of  such  theories  are  cited 
in  [2-5].  The  simplest  of  all  the  improved  laminated  theories  are 
the  ones  based  on  an  assumed  displacement  field  of  the  form 

u - u°  (x,y, t)  ♦ Z'»'x(x,y,t) 

v - v°  (x,y , t)  + Z<Jiy(x,y,t)  (1) 

w - w (x,y,t) 
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Actually,  relations  (1)  apply  to  both  the  classical  form  and  the 

improved  shear  deformation  theory,  the  difference  being  that  in  the 

classical  theory  the  and  terms  are  directly  specified  in  terms 

x y 

of  derivatives  of  w.  Despite  the  increased  generality  of  the  shear 
deformation  theory,  the  related  flexural  stress  distributions  show 
little  improvement  over  those  of  the  classical  laminated  plate  theory. 
It  is  apparent  that  higher  order  terms  are  needed  in  the  power  series 
expansions  of  the  assumed  displacement  field  to  properly  model  the 
behavior  of  the  laminates. 

In  this  paper,  a consistent  high  order  laminated  plate  theory  Is 
derived  for  the  flexural  behavior  of  laminated  plates.  The  following 
displacement  field  is  assumed: 

u - u°  (x.y)  + (x.y)  + z2  Cx  (x,y)  * z3$x  (x,y) 

V ■ V°  (x.y)  + zi|/y  (x.y)  + z2  Cy  (x.y)  ♦ z3$y  (x.y)  (2) 

2 

W - w°  (x,y)  + zip2  (x.y)  ♦ z ^ (x.y) 

The  level  of  truncation  in  Equations  (2)  is  consistent  in  the  sense 
that  the  transverse  shear  strains  due  to  in-plane  displacements  u and 
v are  of  the  same  order  in  z as  those  determined  by  the  transverse 

displacement  w.  This  high  order  laminated  plate  theory  is  an  extension 

of  the  thick  plate  theory  developed  earlier  for  homogeneous  Isotropic 
plates  [I].  The  accuracy  of  the  theory  is  assessed  through  its 
application  to  the  problems  of  a bi-directional  and  an  angle-ply 
laminate  subjected  to  sinusoidal  surface  loading. 
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In  addition  to  surface  loading  problems,  there  are  other  classes 
of  problems  of  practical  interest  which  require  the  use  of  a theory  of 
the  order  of  the  one  given  herein.  For  example,  stress  risers  such  as 
cut-outs,  loaded  holes  and  subsurface  cracks,  and  problems  involving 
the  impact  of  laminates  by  foreign  objects  cause  stress  distributions 
and  localized  stress  gradients  through  the  thickness  of  the  plate  that 
are  essentially  three-dimensional  in  nature  and  require  the  application 
of  a high  order  theory. 

LAMINATED  PLATE  THEORY 

The  theory  used  in  this  paper  is  based  on  the  assumed  displacement 
field  as  given  in  Equations  (2).  The  constitutive  relations  for  any 
layer  of  the  laminate  are  of  the  form 
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where  C.^.  are  the  components  of  the  anisotropic  stiffness  matrix.  It 
should  be  noted  that  all  the  six  components  of  stress  and  strain  tensors 
are  included  in  this  theory. 
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The  governing  equations  pertinent  to  this  theory  are  derived 
using  the  principle  of  stationary  potential  energy  in  the  same 
manner  as  that  in  Ref.  [I].  Eleven  equilibrium  equations  are  obtained 
for  the  determination  of  the  eleven  generalized  displacement  coefficients 
in  Eqs.  (2).  The  details  of  the  derivations  are  omitted  here;  suffice 
it  to  say  in  full  form  the  derivation  is  rather  lengthy.  These 
governing  equations  are  recorded  here,  in  an  archive  journal,  for 
possible  future  use  in  related  problems. 

Governing  equilibrium  equations  of  higher  order  plate  theory  are 
given  by 
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are  the  force  resultants  appropr iate  to  this  theory  and 
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where  n and  t are  the  directions  normal  and  tangential  to  the  edge  of 
the  plate. 

The  equations  and  boundary  conditions  given  in  (4)-(7)  are 
independent  of  the  properties  of  the  materials  of  the  plate  and  hence 
hold  true  for  homogeneous  Isotropic  as  well  as  laminated  plates. 

Mext,  the  governing  equilibrium  equations  are  expressed  in  terms 
of  the  displacement  coefficients. 


Using  Eqns.  (2),  (4),  (5)  and  strain-displacement  relations,  the 
governing  equilibrium  equations  can  be  written  in  the  following  form: 
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where  the  notation 


[ ]. 


xx 


etc.  refers  to  derivatives  of  the  column  matrix 


of  generalized  displacements  shown  In  (8). 

The  coefficients  Aj ^ , B ^ ^ , D.j,  F.j,  H.j,  and  £ are  defined 
as  follows 


(A.jt  B j j , D.j,  j » 8jj»  K. j , j j) 
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'-h/2  (10) 

For  symmetric  laminates,  B.j,  F.j,  and  K.j  are  identically  equal  to  zero. 
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Problems  involving  angle-ply  and  cross-ply  laminates  subjected 
to  static  sinusoidal  loadings  are  next  used  to  assess  the  degree  of 
accuracy  of  this  high  order  laminated  plate  theory. 

NUMERICAL  EXAMPLES 
A.  Angle-Ply  Laminate 

Consider  an  infinite  angle-ply  laminate  of  thickness,  h,  subjected 
to  a pressure  on  the  top  surface  z ■ h/2  of  the  form 

q " qo  sln  T 00 

with  all  other  surface  tractions  identically  equal  to  zero.  The  problem 
can  be  solved  by  assuming  a solution  of  the  form 

u - (aQ  + zaj  + z2a2  + z^a^)  cos 

v - (bQ  + zbj  + z2b2  + z^bj)  cos  ^ (12) 

w ■ (cq  + zcj  + z2c2)  sin 

where  the  constants  a.,  b.,  and  C|  are  determined  by  the  satisfaction 

of  the  governing  equations  given  in  the  Appendix. 

Numberical  results  for  a three-layer  symmetric  laminate  are  shown 
in  Figures  1-A.  The  ply  orientations  and  thickness  are  (+30°,  -30°,  +30’) 
and  (h/4 » h/2,  h/ 4)  respectively.  The  following  properties  are  used  for 
each  ply 


03) 


El  - 25  x 106  psi , Et  - 106  psi 
Glt  - 0.5  x 106  psi,  Gtt  - 0.2  x 106  psi 

VLT  “ VTT  * °‘25 

where  L and  T are  the  directions  parallel  and  normal  to  the  fibers, 
respectively,  and  is  the  Poisson's  ratio  measuring  transverse  strain 
under  normal  stress  parallel  to  the  fibers.  These  are  typical  values 
of  high  modulus  graphite/epoxy  composites.  The  stress  and  displacement 
components  in  Figures  1-4  are  normalized  as  follows 


u 


IOOTu 
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(14) 


for  comparison  with  the  exact  elasticity  solutions  given  by  Pagano  [7]. 

Figures  1 and  2 show  the  flexural  stress  distributions  for  the 
case  L/h  ■ 10  and  l/h  ■ 4,  respectively.  The  agreement  with  exact 
elasticity  solution  is  exceptionally  good  in  the  region  of  high  values 
of  flexural  stresses.  As  the  interface  between  different  layers  is 
approached,  the  stresses  in  the  +30®  is  slightly  different  from  that 
given  by  exact  elasticity  solutions.  However,  such  a slight  difference 
is  immaterial,  especially  at  the  regions  of  low  values  of  stresses. 

Figures  3 and  4 show  the  corresponding  In-plane  displacement  in 
the  x-direction.  As  in  the  case  of  flexural  stresses,  good  agreement 
with  exact  elasticity  solutions  is  observed.  The  solutions  in  Figures  3 
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and  4 reveal  the  necessity  for  modeling  the  nonlinear  distribution 
of  displacements. 

B.  Bidirectional  Laminates 

A more  critical  test  of  the  laminated  plate  theory  can  be  obtained 
by  repeating  the  above  problem  for  a symmetric  bidirectional  laminate. 

In  this  case  a higher  discontinuity  in  material  properties  is  experienced 
at  the  interface  of  different  layers.  Numerical  results  for  a three- 
layered (0°,  90°,  0°)  bidirectional  laminate  are  given  in  Figures  5 and  6 
for  the  flexural  stress  distributions.  The  material  properties  in  each 
layer  are  the  same  as  gWen  in  Eqs.  (13)  and  the  results  are  compared 
with  the  exact  elasticity  solutions  given  by  Pagano  [8]. 

As  in  the  case  of  angle-ply  laminates,  close  agreement  of  the 
numerical  results  with  exact  elasticity  solutions  is  obtained.  The 
relatively  large  discrepancies  in  the  values  of  the  flexural  stresses 
at  the  interface  between  different  layers  is  due  to  the  high  discontinuity 
in  the  values  of  Cjj  across  the  interface  of  different  layers.  However, 
as  before,  such  discrepancies  occur  in  the  regions  of  low  values  of 
flexural  stresses  where  accurate  predictions  of  flexural  stresses  are 
not  important. 

DISCUSSION 

By  comparing  the  results  obtained  with  the  exact  elasticity 
solutions  and  the  classical  laminated  plate  solutions,  it  is  obvious 
that  the  present  high  order  laminated  plate  theory  gives  a much  better 
approximation  to  the  behavior  of  laminated  plates.  This  Is  especially 


true  In  the  case  of  relatively  thick  laminates  where  the  effects  of 
the  transverse  components  of  stress  and  strain  could  not  be  neglected. 

As  in  Ref.  [1],  the  present  high  order  theory  is  expected  to  give 
reasonably  accurate  solutions  for  problems  where  the  characteristic 
length  of  the  loading  pattern  or  the  dimensions  of  the  disturbing 
features,  such  as  cut-outs,  are  of  the  order  of  thickness  of  the  plate. 

The  laminate  stiffness  coefficients  8.^.,  F}  ^ , and  K.  ^ , as  given 
in  Eq.  (14)  represent  the  coupling  between  the  in-plane  and  out-of-plane 
response  of  the  laminate.  In  the  case  of  symmetric  laminates,  it  can 
be  seen  from  Eq.  (14)  that  these  coupling  coefficinets  are  identically 
equal  to  zero.  Thus,  for  symmetric  laminates,  the  governing  eleven 
equilibrium  equations  can  be  separated  into  a set  of  five  second-order 
differential  equations  governing  the  in-plane  displacement  components 
of  the  laminate  as  given  by 

u • u“  (x,y)  + z2Cx  (x,y) 

V * V°  (x,y)  + z2cy  (x.)  (15) 

w - zty*  (x,y) 

and  a set  of  six  second  order  differential  equations  governing  the 
flexural  displacement  components  as  given  by 

u - z<Jix  (x,y)  + z3d>x  (x,y) 

v - zi(»y  (x,y)  + z3$y  (x,y)  (16) 

w - w°  (x,y)  + z2S2  (x,y) 
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For  symmetric  laminates  subjected  to  loadings  which  also  are 
symmetrical  with  respect  to  the  middle  plane  of  the  lamina  % Eqs.  (15) 


provide  an  approximation  to  the  in-plane  response  of  the  laminates, 
with  no  contribution  coming  from  Eqs.  (16). 

Similarly,  if  the  symmetric  laminate  is  subjected  only  to 
anti -symmetri cal  loadings,  Eqs.  (16)  provide  an  approximation  to  the 
flexural  behavior  of  the  laminate.  For  loadings  which  are  neither 
symmetrical  nor  ant i -symmetr ical , terms  from  both  Eqs.  (15)  and  (16) 
will  contribute  to  the  total  response  of  the  laminate.  For  many 
problems,  however,  the  in-plane  contributions  to  the  flexural  behavior 
of  symmetric  laminates  are  small  except  for  very  thick  laminates.  To 
illustrate  this.  Fig.  7 gives  the  flexural  stress  distribution  for  a 
moderately  thick  symmetric  angle-ply  laminate  based  on  Eqs.  (2)  and 
(16)  for  l/h  ■ 4.  As  can  be  seen  from  the  figures,  the  two  different 
stress  distributions  are  extremely  close  to  each  other.  Based  upon 
the  results  obtained  in  Part  I,  it  would  probably  be  at  about  L/h  ■ I 
where  the  contributions  from  both  Eqs.  (15)  and  (16)  are  comparable. 

For  unsymmetric  laminates,  the  coupling  coefficients  B.^,  F„  and 
K„  are  not  identically  equal  to  zero  and,  hence,  the  in-plane  and 
out-of-plane  deformations  are  coupled  with  each  other.  The  behavior 
of  the  laminate  is  thus  determined  by  the  solution  of  a set  of  eleven 
coupled  second  order  partial  differential  equations.  Figure  8 gives 
the  flexural  stress  distributions  for  a [+30°,  -30®,  +30°,  -30°,  +30®, 
-30®]  unsymmetric  laminate  subjected  to  a sinusoidal  loading  as 
specified  by  Eqs.  (11).  The  corresponding  solutions  obtained  by 
using  Eqs.  (16),  i.e.,  neglecting  the  coupling  effects  of  the  in-plane 


displacements  components,  are  also  given.  From  Fig.  8 it  is  seen  that 
the  flexural  stress  distributions  are  strongly  influenced  by  the  effects 
of  in-plane  coupling.  In  fact,  the  use  of  only  the  terms  in  (16)  leads 
to  a solution  which  has  a continuous  flexural  component  of  stress, 
which  as  seen  from  Fig.  8,  is  completely  erroneous.  Thus,  Eqs.  (2) 
must  be  used  to  obtain  accurate  stress  distributions. 

Finally,  it  should  be  mentioned  that  the  present  results  were 
obtained  with  no  recourse  to  shear  correction  factors,  which  are 
commonly  employed  in  laminated  plate  analysis.  As  discussed  in  Chapter  I 
of  this  work,  it  would  be  inconsistent  to  employ  these  factors  with  a 
high  order  theory  of  plate  deformation. 
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FIGURE  CAPTIONS 

Figure  1:  Flexural  stress  distributions  for  [+30°,  -30°]s  angle-ply 
laminate  at  L/h  ■ 10.0. 

Figure  2:  Flexural  stress  distributions  for  [+30°,  -30°]s  angle-ply 
laminate  at  L/h  ■ 4.0. 

Figure  3:  In-plane  displacement  IT  for  [+30°,  "3_0]s  angle-ply 
laminate  at  L/h  ■ 10.0. 

Figure  4:  In-plane  displacement  u for  [+30°,  -30°]s  angle-ply 
laminate  at  L/h  - 4.0. 

Figure  5:  Flexural  stress  distributions  for  [0°,  90°,  0°]  cross-ply 
1 ami  ante  at  L/h  - 10.0. 

Figure  6:  Flexural  stress  distributions  for  [0°,  90° , 0°]  cross-ply 
laminate  at  L/h  * 4.0. 

Figure  7-  Comparison  of  flexural  stress  distributions  due  to 

Eqs.  (2)  and  (9)  for  [+30°,  -30°]s  angle-ply  laminate  at 
L/h  ■ 4. 

Figure  8:  Comparison  of  flexural  stress  distributions  due  to 

Eqs.  (2)  and  (9)  for  unsymmetric  [+30°,  *30°,  +30°,  -30°, 
+30°,  -30°]  angle-ply  laminate  at  L/h  « 10. 
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Chapter  3 

STRESS  SOLUTION  DETERMINATION  FOR 
HIGH  ORDER  PLATE  THEORY 

K.  H.  Lo,  R.  M.  Christensen  and  E.  M.  Wu 

u 


ABSTRACT 


The  high  order  theory  of  plate  deformation  developed  in 
References  [I]  and  [2]  is  further  examined  herein.  Specifically, 
stress  solutions  are  given  and  evaluated  against  exact  elasticity 
solutions  under  stringent  short  wave  length  load  conditions.  By  the 
first  method  the  stresses  are  evaluated  directly  from  the  resulting 
• displacement  solutions.  In  a more  refined  procedure,  the  transverse 

[ 

shear  stresses  and  the  transverse  normal  stress  are  evaluated  by  an 
alternate  equilibrium  method.  The  latter  procedure  is  shown  to  be  more 
accurate  than  the  former. 


INTRODUCTION 

It  has  long  been  recognized  that  classical  plate  theory  must  be 
modified  to  treat  certain  high  order  effects.  The  first  comprehensive 
generalization  of  the  classical  theory  was  that  given  by  Reissner  [3]. 

Since  Reissner* s work,  there  have  been  a great  many  further  generalizations 
beyond  the  classical  theory  assumptions,  with  perhaps  the  highest  order 
theory  to  date  being  that  given  by  Lo,  Christensen  and  Wu  [I]  and  [2]. 
Preliminary  steps  were  taken  in  References  [1]  and  [2]  to  assess  the 
accuracy  of  the  theory.  In  this  paper  this  important  subject  is 
examined  in  greater  detail. 

The  theory  developed  in  References  [1]  and  [2]  is  based  upon  an 
assumed  displacement  field  of  the  type 

u - u®(x,y)  + z^x(x,y)  + z2Cx(x,y)  + z3$x(x,y) 

v - v®(x,y)  + z’^y (x,y)  + z2Cy(x,y)  -t  z3<j>y(x,y)  (1) 

w ■ w® (x,y)  + z’j>2(x,y)  + z2C2(x,y) 
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where  u and  v are  the  in-plane  displacement  components,  w the  out  of 
plane  or  transverse  components,  z the  normal  coordinate,  and  the 
remaining  functions  in  (1)  depend  upon  the  in-plane  coordinate  x and  y. 
The  governing  theory,  based  upon  the  principle  of  stationary  potential 
energy,  resulted  in  eleven  second  order  partial  differential  equations 
to  determine  the  eleven  functions  in  (I).  It  appears  that  an  approach 
of  this  type  is  the  logical  way  to  proceed  if  one  wishes  to  determine 
only  the  displacements.  It  is  less  clear  that  this  approach  is  the 
most  expeditious  method  if  one  seeks  to  determine  stresses.  Now  the 
comparisons  with  exact  elasticity  results  given  in  References  [I]  and 
[2]  were  only  for  the  in-plane  stress  components,  the  transverse  shear 
stresses  and  the  transverse  normal  stress  were  not  evaluated.  Therefore 
the  more  complete  stress  information  to  be  given  here  will  help  to 
answer  the  question  of  the  general  accuracy  of  the  theory.  Before 
proceeding  with  this  however,  it  is  useful  to  consider  the  three 
theoretical  approaches  to  plate  and  shell  development,  and  some 
advantages  and  disadvantages  of  each. 

The  first  and  most  obvious  approach  to  deriving  an  approximate 
plate  theory  utilizes  assumptions  upon  the  forms  of  displacements,  as 
in  (1).  The  governing  differential  equations  could  then  be  derived 
either  by  a direct  method  as  in  the  case  of  classical  plate  theory, 
or  by  the  use  of  the  principle  of  stationary  potential  energy  as  in 
Reference  [1].  Equilibrium  is  violated  by  this  approach,  that  is  to 
say,  the  equilibrium  equations  are  only  approximately  satisfied  through 
weighted  averages.  The  second  possible  approach  is  the  direct  reversal 
of  that  just  described.  Stress  expansions  in  z are  assumed  and  the 
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governing  differential  equations  are  derived  either  by  a direct 
approach  or  by  the  use  of  the  complementary  energy  principle.  Typical 
examples  of  this  approach  include  the  work  of  Reference  [4].  In  this 
method  the  equilibrium  equations  are  satisfied,  the  stress  strain 


relations  are  satisfied,  but  the  compatibility  of  displacements  is 


violated.  In  the  third  approach,  assumptions  are  made  upon  both  the 
stress  states  and  the  displacement  forms  where  by  both  equilibrium 
and  compatibility  conditions  are  satisfied.  However  the  stress  strain 
relations  are  violated.  Reissner's  plate  theory  is  the  most  common 
form  of  this  type,  Reference  [3]. 

It  is  the  stress  state  which  usually  is  the  item  of  interest  in 
most  problems.  Accordingly  it  might  seem  to  be  most  rational  to  use  either 
of  the  latter  two  stress  type  theories,  but  not  the  theory  which  depends 
exclusively  upon  displacement  assumptions.  However,  such  reasoning 
involves  one  serious  problem,  namely  it  would  eliminate  the  extension 
of  the  theory  to  model  laminate  behavior.  This  consequence  is 
because  of  the  fact  that  in  laminates  the  in-plane  stress  components 
are  discontinuous  and  it  becomes  as  exceedingly  complex  matter  to 
construct  a high  order  theory  which  must  inherently  account  for 
discontinuous  stresses.  However,  even  in  laminates  the  displacements, 
of  course,  are  continuous,  and  since  a major  impetus  in  constructing 
new  high  order  theories  is  for  use  with  laminates,  it  is  herein 
considered  necessary  to  proceed  with  the  displacement  theory  of 
Reference  [I].  The  displacement  theory  of  Reference  [I]  in  fact  has 
been  extended  to  model  laminates  by  Lo,  Christensen  and  Wu  [2]. 


! 


I 
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This  status  of  affairs  still  leaves  us  with  some  uncertainty. 

Is  the  displacement  type  theory  of  Referinces  [I]  and  [2]  the  best 
means  by  which  to  deduce  the  stress  state  under  conditions  where 
high  order  effects  are  of  importance?  The  doubt  arises  because  the 
equilibrium  conditions  are  violated  by  the  theory  wheras  stresses 
possess  a one  to  one  correspondence  to  the  equilibrium  conditions. 
However,  there  is  one  possible  means  by  which  the  accuracy  of  the 
stresses  obtained  by  this  displacement  type  theory  can  be  improved. 

The  possible  procedure  is  as  follows.  Use  the  high  order  theory  based 
upon  (1)  to  deduce  the  in-plane  stress  components,  o^,  o^,  and  T^. 

Then  insert  these  stress  solutions  into  the  equilibrium  equations 
and  solve  for  the  out  of  plane/transverse  stress  components  t 
and  oz  by  integration.  This  procedure  clearly  results  in  a stress 
solution  which  satisfies  equilibrium  exactly.  The  procedure  is 
suggested  by  the  classical  theory  approach,  which  does  not  directly 
provide  a solution  for  the  transverse  stress  components,  and  they 
have  to  be  found  by  the  method  described  above. 

Thus,  the  stresses  implicity  in  the  high  order  displacement 
type  theory  of  References  [1]  and  [2]  will  be  determined  by  two 
separate  means.  First  the  in-plane  and  transverse  stresses  will  be 
found  directly  from  the  displacement  solution  through  the  use  of  the 
strain-displacement  and  stress  strain  relations.  By  the  second  method 
the  in-plane  stresses  will  be  found  by  the  method  just  described  and 
they  will  then  be  used  to  determine  the  transverse  stresses  by 
integrating  the  equations  of  equilibrium.  These  two  alternate  methods 
of  deducing  stress  will  be  compared  and  tested  against  exact  elasticity 
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solutions.  Both  cases  of  homogeneous  plates  and  laminates  will  be 
considered. 

THEORETICAL  CONSIDERATIONS 

The  assumed  displacement  fields  are  given  by  relations  (1).  The 
polynomial  expansion  for  w is  truncated  at  one  order  lower  than  the 
expansions  for  u and  v such  that  the  contributions  to  the  transverse 
shear  strains  from  u and  v are  of  the  same  order  in  z as  that  from 
che  terms  in  w.  The  strain-displacement  relations  of  the  linear  theory 
of  elasticity  are 


2 ? 

u® , + ztfi  + z C + zJb 

’x  Yx,x  X,X  YX,X 

V®.  + zl p + z2?  + z^4> 

y y,y  y.y  y»y 


2z'z 


Y 

= Y ® 

+ 

zr 

+ 

z23 

xy 

‘xy 

xy 

xy 

Y 

= Y 8 

+ 

zr 

+ 

z\_ 

'xz 

xz 

xz 

xz 

Y 

= Y ° 

+ 

zr 

+ 

z2B 

yz 

yz 

yz 

yz 

Y • =»  u®,  + v®,  ; y ® " > + w®,  ; y.  + w®, 

Txy  ’y  ’x’  'xz  yx  ’x’  yz  Yy  ’y 

F*y  - *x.,  + *y,x=  rxz  - *z.x=  ryz  - «y  + *z,y 

^xy  " ?x,y  + S,x;  ^xz  = ^x  + ?z,x’  ®yz  = ^y  + ^z,y 

Nxy  x,y  y,x 


« 


The  stress  strain  relations  appropriate  to  an  anisotropic  material 
are  given  by 
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where  C j j , i ,j  » I,  2,  ...,  6 are  the  stiffness  coefficients.  The 
derivation  of  the  governing  equations  for  this  higher  order  plate 
theory  is  given  in  Reference  [1]  for  homogeneous  isotropic  plates, 
and  its  extension  to  laminated  plate  conditions  is  given  in 
Reference  [1*]  . The  evaluation  of  the  stresses  will  now  be  given  by 
the  two  methods  mentioned  In  the  introduction. 


STRESS  EVALUATION.  HOMOGENEOUS  PLATE 

An  infinite  homogeneous  isotropic  plate  is  subjected  to  sinusoidal 
loadings  as  in  Reference  [1].  The  prescribed  surface  tranctions  are 

°z(h/2)  " Vin  (T-  V“h/2)  * 0 (6) 

and 
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Txz  h/2^  " Tyz  h/2^  ■ 0 
where  h is  the  plate  thickness. 
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The  solutions  to  the  plate  problems  are  [I], 


u " I*,  + *A,  + z2A2  + z3A3]  cos  (^) 


V * O 


W - [Cj  + 2C^  + 2 C^]  sin  (~) 

where  the  constants  Aq,  Aj , A2>  A^  and  Cj,  C2,  and  Cj  are  obtained  by 
satisfaction  of  the  governing  differential  equations  and  boundary 
conditions.  From  (3),  (b)  and  (5)  the  transverse  stress  components 
appropriate  to  a homogeneous  isotropic  plate  are 
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Substituting  (8)  into  (9)  and  (10)  gives  the  transverse  stress 
components  for  the  infinite  plate 
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These  stresses  are  of  course  those  deduced  directly  from  the  theory,  and 
they  will  be  displayed  graphically  for  particular  values  of  h/L,  i.e. 
the  ratio  of  the  thickness  of  the  plate  to  the  characteristic  length  of 
the  loading  pattern.  First  of  all,  however,  these  same  stress  components 
will  be  evaluated  by  the  alternate  method  mentioned  in  the  introduction. 

The  in-plane  stress  components  for  the  problem  under  consideration 
are  given  in  Reference  [1],  These  in-plane  stresses  are  substituted  into 
the  equilibrium  equations 
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:o  yield,  after  integration,  the  out  of  plane/transverse  stress 
:omponents: 
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where  the  top  and  bottom  surface  traction  conditions  (6)  and  (7)  have 
been  used  to  evaluate  the  constants  of  integration. 

The  corresponding  exact  elasticity  solution  for  this  problem  is 
given  in  Reference  [5]*  among  many  other  sources.  Before  comparing 
the  two  alternate  means  of  deriving  the  stresses,  first  the  corresponding 
results  for  a particular  laminated  plate  will  be  stated. 

LAMINATED  PLATE 

The  laminate  to  be  considered  is  that  of  symmetric  cross-ply 
geometry  with  each  lamina  being  orthotropic.  Using  relations  (5), 

(8),  and  (13),  the  transverse  stress  components  are  found  to  be 


where  the  index  k refers  to  the  k1^1  layer  and  where  f (x)  and  g^(x) 

are  determined  from  the  boundary  and  continuity  conditions.  These 

transverse  stresses  are  those  found  by  the  method  involved  in  using  the 

in-plane  stresses  in  the  equations  of  equilibrium  which  are  then 

integrated  to  find  the  three  remaining  transverse  stress  components. 

Further  consideration  will  be  limited  to  the  case  of  a three  layer 

laminate  arranged  with  the  fiber  directions  designated  by  [0°/90o/0°] 

and  with  the  laminae  having  equal  thickness  of  h/3.  The  following 

fkl  k 

expressions  are  obtained  for  f'  (x)  and  g (x) . 
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where  the  subscripts  I,  2,  and  3 refer  to  the  0°,  90°,  and  0°  layers, 
respectively. 

The  stress  distributions  will  be  displayed  for  lamina  of  the 
following  properties  which  are  typical  of  high  modulus  graphite/epoxy 


compos i tes 


E2  - 25  x I06  psi;  Et  - I06  ps i 
GlT  - 0.5  x 106  psi;  G.^  - 0.2  x 106  psi 


VLT  " VTT  “ °'25 

where  L and  T refer  to  the  properties  alone  and  transverse  to  the  fiber 
directions,  respectively,  and  v^.  ?s  the  Poisson's  ratio  measuring 
transverse  strain  under  normal  stress  parallel  to  the  fibers. 


DISCUSSION 

The  transverse  normal  stress  and  transverse  shear  stress  for  the 
homogeneous  isotropic  plate  case  are  given  in  Figures  1-3.  The 
difference  in  these  figures  is  due  to  the  variation  in  the  ratio  of 
h/L,  i.e.  the  ratio  of  thickness  to  half  wave  length  of  the  sinusoidal 
load.  In  Figure  1,  for  h/L  ■ 1/1*,  the  normal  stresses  calculated  by 
the  two  different  means  are  compared  with  the  exact  solutions.  For 
this  small  value  of  h/L,  the  transverse  shear  stresses  calculated  by 
the  two  methods  are  so  close  to  the  exact  solutions  as  to  make  them 
look  identically  equal  graphically.  Clearly  the  transverse  normal 
stresses  found  by  the  integration  of  the  equations  of  equilibrium  are 
far  more  accurate  than  those  obtained  directly  from  the  displacement 
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solution  through  the  strain  displacement  and  the  stress  strain  relations. 
The  results  Shown  in  these  figures  reveal  that  the  transverse  stresses 
found  by  equilibrium  equation  integration  to  be  more  accurate  than  those 
found  directly  from  the  displacement  solution,  with  these  results  being 
under  stringent  short  wave  length  load  conditions.  As  discussed  in 
Reference  [1]  the  maximum  ratio  of  h/l  for  which  the  theory  has 
reasonable  validity  is  about  h/L  ■ I and  the  results  shown  here 
collaborate  this  conclusion. 

It  is  of  interest  to  note  from  Figures  1-3  that  the  transverse 
stresses  obtained  directly  from  the  displacement  solution  violate  the 
top  and  bottom  traction  conditions.  An  examination  of  the  derivation 
in  Reference  [1]  reveals  this  to  be  a consistent  aspect  of  the  method. 
Thus  even  though  the  tractions  enter  the  theory  as  boundary  conditions, 
this  process  actually  occurs  through  an  equilibrium  weighting  method, 
thus  the  theory  does  not  provide  exact  satisfaction  of  these  boundary 
conditions.  Consider  however,  the  alternate  method  of  obtaining  the 
transverse  stresses  from  integrating  the  equilibrium  equations 
utilizing  the  in-plane  stresses  found  directly  from  the  displacement 
solution.  In  this  case  the  boundary  tractions  are  automatically 
satisfied  through  the  evaluation  of  the  constants  of  integration.  A 
similar  situation  exists  in  the  case  of  laminated  plates.  Transverse 
stresses  evaluated  directly  from  the  displacement  solution  would  in 
general  not  be  continuous  across  the  interfaces  between  lamina, 
however  the  transverse  stresses  found  by  the  equilibrium  method 
proposed  herein  provide  continuous  stress  with  exact  satisfaction  of 
top  and  bottom  surface  conditions. 
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The  results  for  a three  layer  laminate  are  shown  in  Figures  k 
and  5.  for  the  value  of  h/L  ■ 1A.  The  exact  elasticity  solution  is 
taken  from  Pagano  [6].  It  is  apparent  that  the  case  of  the  laminate 
provides  much  more  strenuous  conditions  against  which  to  test  a plate 
theory  than  does  homogeneous  conditions.  Nevertheless,  as  seen  in 
these  figures  the  equilibrium  equation  method  of  generating 
transverse  stresses  provides  a reasonable  approximation  to  the  exact 
solution. 

CONCLUSIONS 

The  present  high  order  theory  of  plate  deformation  appears  to 
provide  reasonably  accurate  predictions  of  behavior  under  short  wave 
length  conditions.  This  conclusion  is  valid  for  both  homogeneous 
plates  and  for  laminates;  also  as  shown  by  the  results,  laminates  are 
much  more  demanding  of  high  order  effect  representation  than  are 
homogeneous  plates.  In  problems  where  displacements  are  the  quantity 
of  prime  interest  the  present  displacement  type  theory  appears  to  provide 
a reasonable  and  high  order  effect  solution.  In  problems  where  the 
stresses  are  the  quantity  sought  it  has  been  shown  that  the  present 

HI 

theory  still  provides  highly  accurate  stress  information.  It  has 
been  demonstrated  that  the  best  method  for  determining  the  stresses 
involves  determining  the  in-plane  stresses  directly  from  the  displacement 

1 

solution  and  thence  determining  the  transverse  stresses  through  the 
integration  of  the  equations  of  equilibrium,  utilizing  the  in-plane 

I 

solution  therein.  This  method  is  of  course  applicable  to  a theory 
of  any  order  not  just  the  present  high  order  theory.  The  success  of 
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the  method  was  virtually  assured  by  the  fact  that  It  is  the  only 
possible  procedure  for  use  at  the  level  of  the  classical  theory 
assumptions . 

I 

1 
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FIGURE  CAPTIONS 


Figure  1. 


Figure  2. 


Figure  3. 


Figure  4. 


Figure  5. 


Transverse  normal  stress  distributions  for  a homogeneous 
isotropic  plate  at  h/L  * 0.25. 

Transverse  normal  stress  distributions  for  a homogeneous 
isotropic  plate  at  h/L  * 1.0. 

Transverse  shear  stress  distributions  for  a homogeneous 
isotropic  plate  at  h/L  * 1.0. 

Transverse  normal  stress  distributions  for  a [0°,  90°,  0°] 
laminate  at  h/L  = 0.25. 

Transverse  shear  stress  distributions  for  a [0°,  90° , 0°] 
laminate  at  h/L  * 0.25. 
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F i gure 


Figure  *». 


Chapter  k 


THREE  DIMENSIONAL  FAILURE  CRITERION 


E.  M.  Wu,  R.  M.  Christensen  and  K.  H.  Lo 


I 


f 


Strength  of  Thick  Laminates 

The  strength  of  thick  laminates  can  be  cast  in  terms  of  the 
characterization  of  the  criticality  of  three-dimensional  stress  states, 
which  can  lead  to  rupture.  Due  to  the  complexity  of  this  characterization, 
the  following  simplified  assumptions  are  required:  1)  all  combined 
loadings  are  proportional  and  monotonic  thus  precluding  the  effects  of 
fatigue  and  loading  path  dependency,  2)  the  thick  composite  is  planar 

I* 

laminated  thus  precluding  3-D  woven  structures.  Within  these  constraints, 
significant  flexibility  is  retained,  e.g.  the  lamination  sequence  and 
angle  can  be  arbitrarily  varied  utilizing  a single  set  of  basic  data, 
or  if  basic  data  are  available  for  several  generlcally  different  laminae, 
thick  hybrid  composite  laminate  strengths  can  also  be  estimated. 

The  flexibility  of  this  characterization  lies  in  the  inclusion  of 
three-dimensional  stresses  in  the  tensor  polynomial  failure  criterion 
for  the  lamina.  Since  the  coefflcinets  of  the  tensor  polynomial 
failure  criterion  obey  tensor  transformation  laws,  computational  methods 
for  predicting  the  strength  of  thick  laminae  of  any  layup  can  be 
establ i shed. 

The  principle  steps  of  this  characterization  are: 

1)  Evaluation  of  the  3“D  failure  criterion  for  a unidirectional 
lami na, 

2)  Establishment  of  the  transformational  properties  of  the  strength 
coefficients, 

3)  3"0  analysis  of  thick  laminated  structure. 

These  steps  are  illustrated  in  Figures  la,  b,  c respectively.  The  last 
step  has  already  been  addressed  in  Chapters  I,  2,  and  3.  The  first  and 
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second  steps  will  now  be  examined.  In  Step  1 (Figure  la)  we  seek  the 
strength  response  of  the  lamina  (along  its  principal  directions)  which 
can  be  expressed  in  terms  of  the  tensor  failure  polynomial  criterion  as  [I]: 
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i.j  - 1 ,2,3, 4,5.6 
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The  failure  tensors  which  characterize  the  strength  of  the  laminate  can 
be  written  in  matrix  form,  while  accounting  for  orthotropy,  as 
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It  is  seen  from  Equation  (lb)  that  the  failure  tensors  associated  with 
the  3 axis  (F^,  F^,  F^,  F^,  F^,  F^)  control  the  strength  related 
to  the  stresses  assocated  with  the  3 axis.  These  coefficients  are  not 
required  if  the  laminate  is  thin  but  must  be  included  if  the  laminate 
is  thick.  The  remaining  coefficlnets  are  well  known  and  they  may  be 
measured  by  a total  of  six  independent  experiments.  It  appears  that 
measurement  of  the  failure  tensors,  associated  with  the  inclusion  of  the 
3 axis  requires  six  additional  tests.  However  symmetry  conditions  will 
be  explored  to  reduce  the  number  of  tests. 

A unidirectional  lamina  (oriented  as  shown  in  Figure  1(a))  possesses 
the  symmetry  condition  that  properties  associated  with  the  are 


82 


V . 

interchangeable  with  properties  associated  with  the  axis.  Thus,  the 
following  table  of  symmetry  for  the  strength  coefficients  can  be  constructed. 

Table  I 


Tensor  Notation 

Contracted  Notation 

F22  " F33 

F2  = F3 

F 1 1 22  = Fl  1 33 

F12  " F1 3 

F a F 

2222  3333 

F22  = F33 

F a F 

1313  1212 

F55  = F66 

By  virtue  of  these  symmetries,  we  can  assess  the  necessary  experiments 
required  to  evaluate  the  coefficinets  for  a 3-0  failure  criterion. 

These  experiments  can  be  determined  by  evaluation  of  Equation  (la) 
following  the  procedure  outlined  by  Wu  in  Reference  [2].  By  substitution 
of  different  stress  states  in  Equation  (la),  we  can  evaulate  the 
coefficients  of  the  strength  tensors: 


Evaluation  of  Component  State  of  Stress 

of  Strength  Tensor* Requ  i red 


\ 

( 


Here,  we  adopt  the  same  notation  as  in  Reference  [2]  where  X. 
represents  the  uniaxial  strength  induced  by  the  uniaxial  stress  o. 
and  o.  are  the  biaxial  strengths  induced  by  combined  stress  conditions. 
The  strength  tensor  is  computatable  from  the  strengths  by  solution 
of  Eq.  (1).  The  necessary  relations  are:  for  experiments  No.  1,  2,  3, 
4,  5,  and  8, 


I,  2,  3,  4,  5,  6. 


for  experiments  No.  6, 


F!2  ' " - F131  - F232  - Fll5!2  ' F22322) 

20,  a2 


For  experiment  No.  7, 


F23  " ' F2a2  " F3°3  " F22J2  * F33°3  ) - (i,) 

We  note  that  tests  No.  (1)  through  (6)  are  those  required  for  a 
two-dimensional  failure  criterion;  the  procedure  for  determining 
them  has  been  thoroughly  described  in  Reference  [2].  In  addition,  we 
uncovered  a rather  surprising  and  convenient  fact  that  a three- 
dimensional  failure  criterion  (in  2nd-order  form)  only  required  two 
additional  experiments  i.e.,  tests  No.  7 and  No.  8.  In  fact,  we  note 
that  experiemnts  No.  7 and  No.  8 are  measuring  the  strength  response  of 


85 


essentially  isotropic  properties.  Thus,  a not  unreasonable  estimate 
of  these  strengths  may  be  arrived  at  by  degrading  the  correspond ing 
resin  matrix  strengths  by  a suitable  factor  to  reflect  the  stress 
concentration  due  to  the  presence  of  fibers.  If  such  a comprise  is 
deemed  acceptable;  ^additional  tests  are  required  to  extend  a 
2-0  failure  criterion  to  a 3“D  failure  criterion!  Such  is  the 
remuneration  for  the  tensor  polynomial  failure  criterion  formulation. 
In  this  program  these  coefficients  are  experimentally  measured  rather 
than  estimated. 


Experimental  Measurements 

The  failure  criteria  outlined  here  are  evaluated  for  a 
ultra-high  stiffness  graphite  epoxy  compos i te-GY  70  manufactured  by 
Celanese  Corporation  with  Fiberite  931*  Resin.  Samples  for  this 
program  was  fabricated  by  General  Dynamics  Cooperation,  Convair 
Division. 

The  eight  characteristic  experiments  described  previously  were 
carried  out.  The  sample  configurations  of  these  tests  are  as  follows. 

Experiment  No.  1 (longitudinal  tensile  test)  was  performed  on 
para  1 lei -edged  samples  with  end-tabs  of  the  configuration  referred 
to  as  IITRI  samples  recommended  by  the  Air  Force  Design  Guide. 

End  tabs  were  not  tapered  and  were  fully  loaded  by  wedge  grips  over 
its  entire  length  to  eliminate  peeling  by  tensile  stress. 

Experiment  No.  2 (longitudinal  compression)  was  performed  on  a 
sandwich  specimen.  The  configuration  is  shown  in  Fig.  2a  and  b. 


This  compression  sample  is  similar  to  the  honeycomb  sandwich  sample 
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suggested  in  the  Air  Force  Design  Guide  (M I L- STD-40 > A)  with  the 
difference  that  low  density  rigid  foam  is  used  in  lieu  of  the  usual 
honeycomb  for  expediency  of  sample  preparation.  Test  results  are 
consistent  with  those  produced  by  the  honeycomb  sandwich  samples. 

Experiment  No.  3 (Transverse  tension)  samples  were  identical 
to  experiment  No.  I with  the  exception  that  end  tabs  were  not  used. 

Experiment  No.  A (Transverse  compression)  samples  were  identical 
to  those  used  in  experiment  no.  2. 

Experiment  No.  5 (Longitudinal  shear)  was  performed  by  a beam 
bending  configuration  as  shown  in  Fig.  3a  and  3b.  This  shear 
test  was  originally  suggested  by  Messmer  [3]  in  photoelasticity 
experiments  and  subsequently  rediscovered  by  losipescu  [4]  to  test 
shear  strength  of  isotropic  metals.  The  salient  features  of  this 
configuration  are  that: 

i)  The  one  dimensional  nature  of  the  configuration  lends  itself 
to  analysis  which  is  accurate  to  the  beam  theory  level. 

ii)  In  beam  bending  theory,  longitudinal  shear  is  proportional 
to  the  gradient  of  the  bending  moment;  in  this  configuration  the  bending 
moment  gradient  is  uniform  over  the  length  of  the  beam  and  the 
magnitude  of  the  moment  is  zero  at  the  center.  Thus  if  failure  is 
induced  to  occur  at  the  center,  the  strength  measurement  reflects  the 
pure  shear  strength.  For  the  purpose  of  inducing  shear  failure  at 
the  center  and  for  the  purpose  of  converting  the  parabolic  distribution 
of  shear  stress  to  uniform  shear  stress,  Losipesu  [4]  introduced  90° 
notches  to  the  sample.  While  this  technique  was  reported  to  work 
satisfactorily  for  metals,  it  proved  to  be  unsuccessful  for 
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unidirectional  lamina  because  machining  of  the  90°  notch  inevitably 
induced  cracking.  We  overcame  this  difficulty  by  introducing 
bonding  flanges  to  the  composite  web  to  form  the  shape  of  an  I beam. 

The  flanges  (top  and  bottom)  are  discontinuous  at  the  center  (see 
Figure  4);  since  the  moment  is  zero  in  the  center  with  no  normal  stress 
due  to  bending,  no  stress  singularity  is  induced  by  the  slit.  At  the 
same  time,  shear  failure  can  propagate  through  the  slit.  The  absence 
of  singularity  is  valid  to  the  level  of  the  strength  of  material 
analysis  but  probably  not  in  the  theory-of-elast ici ty  level.  However, 
verification  was  carried  out  to  the  extent  that  literature  values  on 
known  composites  (Scotch-ply-1002)  were  recovered  by  this  test. 

It  appears  this  test  might  merit  more  extensive  development  under  a 
different  scope.  Suggested  studies  are  finite  element  and  photo- 
elastic analysis  and  optimization  of  the  flange  material. 

Test  No.  6 involved  internal  pressurization  and  compression  of  a 
thin-walled  tube.  The  stress  ratio  induced  was  B ■ o^/c^  * 

Experimental  details  are  as  described  in  Reference  (5). 

Test  No.  7 was  similar  to  a constrained  compression  test 
reported  by  T.  C.  Collings  [5]  and  is  schematically  illustrated  in 
Figure  5*  The  sample  is  a 3/4- inch  cube  cut  from  a thick  lamina.  For 
the  measured  Poisson  ratio,  ■ 0.62,  the  stress  ratios  are 

°_2  -1 
o3  " hT?2  ' 

Test  No.  8 was  performed  in  the  same  fixture  as  in  Test  No.  5 
(the  longitudinal  shear  test),  the  exception  being  that  in  this  case  the 
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fiber  orientation  is  perpendicular  to  the  web  (See  Figure  6).  These 
samples  are  sliced  from  a 0.75"-thick  laminate  with  a diamond  circular 
cutter. 

All  of  the  above  experiments  were  performed  on  an  electro-hydraulic 
servo-controlled  testing  machine.  Load  is  used  as  the  feedback  control 
and  the  loading  function  is  constant  loading  rate  to  failure.  Loading 
rates  for  each  type  of  tests  were  adjusted  (by  trial)  such  that  starting 
from  no  load,  the  time  to  fail  was  approximately  10  minutes.  The 
experimental  results  are: 


Test 

1.  Longitudinal  Tension 

2.  Longitudinal  Compression 

3.  Transverse  Tension 

A.  Transverse  Compression 

5.  Longitudinal  Shear 

6.  Biaxial  Strength  1,  2 plane 

7.  Biaxial  Strength  2,  3 plane 


Transverse  Shear 


Strength  (ksi) 


X' 

A1 


X' 

*2 


(1st  Batch) 

(2nd  Batch) 


Range 

55-5  - 101 
60.9  - 98.6 
69-3  - 91.2 
A. 2 - 5.6 
8.1  - 9.0 
7.3  - 8.4 
7.9  - 8.5 


8. 


Two  batches  of  samples  were  tested  to  assure  consistency  of  material 
and  instrumentation  reproducibility.  Stress-strain  relations  were 
recorded  and  are  shown  in  Figures  7 to  1 1 . In  Figure  7,  the  individual 
stress-strain  curves  for  longitudinal  tension  are  presented.  In 
Figure  8,  these  tests  are  statistical ly  averaged;  the  curves  shown 
are  the  averaged  stress-strain  curve  and  the  + Due  Standard  deviation 
of  the  stress-strain  curves;  "X"  are  the  ultimate  values  for  each 
individual  samples.  In  similar  forms.  Figures  9 and  10  depict  results 
of  transverse  compression  tests.  Figures  11a  and  lib  are  respectively 
the  stress-strain  components  in  the  2 and  3 directions  for  biaxial 
comoression. 

Based  on  these  measured  strengths  the  strength  tensor  is  computed 
in  accordance  with  equations  2,  3 and  4;  together  with  the  symmetry 
conditions,  we  deduced  the  strength  tensor  for  GY  70/934  to  be: 


F;  - 0.157  x 10  J 

F2  - 74.8  x 10"3 

F,  - (74.8  x 10'3) 


Fn  - 0.150  x 10  J 
F22  - 23.8  x 10'3 


(23.8  x 10'J)  I , 
0.954  x 10'3  f ks i 2 


F)3  - (0.954  x 10"3) 
F,,  - -0.027 


The  bracketed  values  are  inferred  from  symmetry  conditions. 


Discussion  and  Conclusion 

With  these  strength  tensor  components,  we  have  a full  characterization 
of  the  3-0  strength  of  the  GY  70/934  composite.  This  3"D  failure  criterion 
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may  be  used  to  assess  the  criticality  of  a three-dimensional  state  of 


stress  either  computed  from  3"dimensional  analysis  or  the  aforepresented 


high-order  plate  analysis.  Furthermore  we  note  that  biaxial  compressive 


stress  tends  to  increase  the  strength  of  the  composite.  For  example  both 


the  biaxial  tests  No.  6 and  7 produced  strength  greater  than  the  uniaxial 


strength  and  respectively.  This  is  reflected  in  positiveness  of 
the  coefficient  of  strength  tensor  F^>  F23*  Physically  this 


suggests  that  the  composite  test  is  probably  microflaw  sensitive  in  the 


transverse  direction;  application  of  transverse  compression  tends  to 


retard  flaw  growth  and  lead  to  a strength  increase  in  the  longitudinal 


direction.  This  may  be  the  reason  for  the  large  range  in  strength 


variability  measured.  Constraints  in  the  program  preclude  further 


in-depth  experimental  measurement.  Another  effort  has  been  initiated  to 


check  the  experimental  technique  and  provide  further  measurements. 


Nevertheless,  the  methodology  and  analysis  remain  well  founded  and  may 


be  applied  with  confidence. 
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1:  Composite  Laminates 

2:  Compression  Sample  Configuration  (Test  No.  2 and  A) 

3:  Shear  Fixture  Configuration  (Test  No.  5) 

A:  Shear  Sample  Configuration  (Test  No.  5) 

5:  Biaxial  Compression  Configuration  (Test  No.  7) 

6:  Shear  Test  Configuration  (Test  No.  8) 

7:  Individual  Stress-Strain  Curves  of  Longitudinal  Tension 

Tests  of  GY-70/93A  (I  ksi  - 0.IA5  MPa)' 

8:  Averaged  Stress-Strain  Curves  of  Longitudinal  Tests  of 

GY-70/93A  (I  ksi  - 0.145  MPa) 

9:  Individual  Stress-Strain  Curves  for  Transverse  Compression 

Tests  of  GY-70/93A  (I  ksi  - 0.145  MPa) 

10:  Averaqe  Stress-Strain  Curves  of  Transverse  Compression 

Tests  of  GY-70/93A  (I  ksi  - 0.145  MPa) 

I la:  Individual  o2  vs  £j  Curves  for  Biaxial  Compression  in 
2-3  plane  for  GY-70/93A  (I  ksi  - 0.145  MPa) 

lib:  Individual  Oj  vs  £j  Curves  for  Biaxial  Compression  in 

2-3  Plane  for  GY-70/93A  (1  ksi  - 0.145  MPa) 
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Figure  2.  Compression  Sample  Configuration  (Test  No.  2 and  4) 
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Figure  7 


Individual  Stress-Strain  Curves  of  Longitudinal  Tension  Tests  of 
GY-70/931*  (1  ksi  - 0.U5  MPa). 
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Averaged  Stress-Strain  Curves  of  Longitudinal  Tests  of  GY-TO/^** 
(1  ks i - 0.145  MPa). 
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Figure  9 


Individual  Stress-Strain  Curves  for  Transverse  Compression  Tests 
of  GY-70/934  (1  ksi  - 0.145  MPa) . 
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Figure  10.  Average  Stress-Strain  Curves  of  Transverse  Compression  Tests 
of  GY-70/93i*  (1  ksi  - 0.1A5  MPa). 
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Figure  I la.  Individual  o2  vs  Curves  for  Biaxial  Compression  In  2-3  Plane 
for  GY-70/931*  0 ksi  - 0.145  MPa). 
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Figure  lib.  Individual  n.  vs  Curves  for  Biaxial  Compression  In  2-3  Plane 
for  GY-70/931*  (I  ksl  - 0.145  MPa). 
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